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Abstract. A new stochastic cellular automaton (CA) model of traffic flow, 
which includes slow-to-start effects and a driver's perspective, is proposed by 
extending the Burgers CA and the Nagel-Schreckenberg CA model. The flow- 
density relation of this model shows multiple metastable branches near the 
transition density from free to congested traffic, which form a wide scattering area 
in the fundamental diagram. The stability of these branches and their velocity 
distributions are explicitly studied by numerical simulations. 



1. Introduction 

Traffic problems have been attracting not only engineers but also physicists 
Especially it has been widely accepted that the phase transition from free to congested 
traffic flow can be understood using methods from statistical physics |21 El ■ In order 
to study the transition in detail, we need a realistic model of traffic flow which should 
be minimal to clarify the underlying mechanisms. In recent years cellular automata 
(CA) 01 have been used extensively to study traffic flow in this context. Due to 
their simplicity, CA models have also been applied by engineers, e.g. for the simulation 
of complex traffic systems with junctions and traffic signals [S]. 

Many traffic CA models have been proposed so far 12 El |S] , and among these 
CA, the deterministic rule-184 CA model (R184), which is one of the elementary CA 
classified by Wolfram is the prototype of all traffic CA models. R184 is known to 
represent the minimum movement of vehicles in one lane and shows a simple phase 
transition from free to congested state of traffic flow. In a previous paper (HI, using the 
ultra-discrete method TIT , the Burgers CA (BCA) has been derived from the Burgers 
equation 

Vt = 2vVx + Vxx, (1) 

which was interpreted as a macroscopic traffic model The BCA is written using 
the minimum function min by 

C/^+i = [/j + min{C/j_i, L ~ Uj} - min{C/j, L - f/j+J, (2) 
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where C/* denotes the number of vehicles at the site j and time t. If we put the 
restriction L = 1, it can be easily shown that the BCA is equivalent to R184. Thus 
we have clarified the connection between the Burgers equation and R184, which offers 
better understanding of the relation between macroscopic and microscopic models. 

The BCA given above is considered as the Euler representation of traffic flow. As 
in hydrodynamics there is an another representation, called Lagrange representation 
|12| . which is specifically used for car- following models. The Lagrange version of the 
BCA is given by [T3| 

xl+^ =xl+mm{Vmax,xl+s-^l-S}, (3) 
where Vmax = S = L and x* is the position of i-th car at time t. Note that in © 
S corresponds a "perspective" or anticipation parameter ^1] which represents the 
number of cars that a driver sees in front, and Vmax is the maximum velocity of 
cars. Q is derived from the BCA mathematically by using an Euler-Lagrange (EL) 
transformation which is a discrete version of the well-known EL transformation 
in hydrodynamics. 

In this paper, we will develop the BCA to a more realistic model by 
introducing slow-to-start (s2s) effects ^1 El ^1 E| and a driver's perspective S. 
Moreover, a stochastic generalization is also considered by combining it with the 
Nagel-Schreckenberg (NS) model piT^. 

2. Traffic models in Lagrange form 

First, let us extend ||2J) to the case Vmax ^ S and combine it with the s2s model. The 
s2s model jj^l is written in Lagrange form as 

=x\+ min{l, - X* - I, x\-\ - x\-^ - I}. (4) 
Note that the inertia effect of cars is taken into account in this model. Comparing 
and (O, we see that, in the s2s model, the velocity of a car depends not only on the 
present headway d* = x-^^— x* — 1, but also on the past headway = x*^J— a;*~^ — I. 
This rule has only a nontrivial effect if d*"^ = and d* = 1, i.e. if the leading car has 
started to move in the previous time step. In this case the following car is not allowed 
to move immediately (s2s). 

Before combining (0) and Q), it is worth pointing out that we can choose the 
perspective parameter as S* = 2 in the model according to observed data. We define 
the size of a cell as 7.5 m and Vmax = 5 in our model according to the NS model. 
Since Vmax corresponds to about 100 km/hour in reality, then one time step in the 
CA model becomes 1.3 s. Moreover, the gradient of the free line and jamming line in 
the fundamental diagram, which is the dependence of the traffic flow Q on density p, 
is known to be about 100 km/hour and —15 km/hour [3] according to many observed 
data (see Fig. Pl [201 1^ . These values correspond to the typical free velocity and the 
jam velocity, respectively. Thus, considering the fact that the positive and negative 
gradient of each line is given by Vmax = 5 and —S/2, respectively, in the CA model 
|12|. we should choose S = 1.5 in the CA model. Since only integer numbers for S 
are allowed in this model, we will simply choose S" = 2 for studying the effect of the 
perspective of drivers. It is noted that other possibilities, such as velocity-dependent 
S or stochastic choice of 5, are also possible. 

Now by combining Q and we propose a new Lagrange model with S ~ 2 
which is defined by the rules listed below: Let vf^ be the velocity of the i-th car at a 
time t. The update procedure from i to t -|- 1 is divided into five stages: 
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Figure 1. An observed fundamental diagram at the Tomei expressway in Japan. 
The gradient of the free Hne (A — B) and jamming Hne (C — D) is known to be 
about 100 km/hour and —15 km/hour. We also see that there is a wide scattering 
area near the phase transition region from free to jamming state. 



1 Accerelation 

vl^^ =mm{Vmax,vf^ + 1}. (5) 

2 Slow-to- accelerate effect 

«f)=min{x;«,xy-xr-2}. (6) 

3 Deceleration due to other vehicles 

z;f)=min{t;f ,a;*+2-^*-2}. (7) 

4 Avoidance of collision 

«f = min{z;f\x*+i-x*-l + z;l^\} (8) 

5 Vehicle movement 

^'=<+^. (9) 

The velocity w^'*^ is used as vf^ in the next time step. ||SJ| is the condition that 
the i-th car does not overtake its preceding (i + l)-th car, including anticipation. The 
accerelation (O is the same as in the NS model, which is needed for a mild accerelating 
behaviour of cars. In the step 2, we call © as "slow-to-accelerate" instead of s2s. This 
is because this rule affects not only the behaviour of standing cars but also that of 
moving cars, which is considered to be a generalization of s2s rule. 

It is not difficult to write down the new model in a single equation for general S. 
The result is 

= X* + min min - x\ - k + V^Vfe)} . (10) 

where the last term represents the collision-free condition explained in Fig. |21 and 

Vl = min{T/„a., - - 5, x\+s -x\-S,x\- x'~^ + 1}. (11) 
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The condition that there is no coUision between the i-th and i + fc-th cars (fc = 
1, • • • , 5* — 1) is given by 

^l+k 'x\~k + y.V, > F/, (12) 

for > 2 (if S" = 1 then we simply put k = V), which is identical to the last term 
in H1U|I . In contrast to the NS model, the velocity of the preceeding car is taken into 
account in the calculation of the safe velocity in the step 4, i.e. our model also includes 
anticipation effects. 



i+k 



i+S 



Vi 



i+k 



Figure 2. Collision-free condition between i-th and i + fc-th car. 



3. Metastable branches and their stabiUty 

Next, we investigate the fundamental diagram of this new hybrid model. In Fig. |31 
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Figure 3. Fundamental diagram of the new Lagrange model. Parameters are set 
to Vmax = 5 and S = 2, and the spatial period is 100 sites. The initial car density 
is varied from 0.05 to 0.95 in steps of 0.01. At each density, we start calculations 
from 30 randomly generated initial configurations, and show only the data at the 
time t = 100. We observe several metastable branches in the deterministic case. 
The fluctuations of the branches show the fact that the asymptotic flow of the 
system sometimes becomes periodic instead of stationary between 0.2 < p < 0.5. 



we observe a complex phase transition from a free to congested state near the critical 
density 0.2 ~ 0.4. There are many metastable branches in the diagram, similar to 
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our previous models in Euler form [221 \2'6\ or in other models with anticipation j24 | . 
We also point out that there is a wide scattering area near the critical density in the 
observed data (Fig. ^ which may be related to these metastable branches. As we 
will discuss later, these branches may account for some aspects of the scattering area 
observed empirically. 

First, we discuss properties of the state in the metastable branches. In all cases it 
consists of pairs of vehicles that move coherently with vanishing headway (see Fig.0J). 
Cars are represented by black squares, and the direction of the road is horizontal right 
and time axis is vertical down. The corresponding velocity distributions are also given 
in Fig.[S| We see that there are stopping cars which velocity are zero only in the case 
of the lowest branch given in the state in Fig. 0] (e) . 




Figure 4. Spatio-temporal patterns of evolution of the uniform flow 
■ ■ ■ 11000001100000 ■ ■ ■ at density p = 2/7 with different strength of perturbation: 
(a) very weak, (b) less weak, (c) medium, (d) stronger and finally (e) strongest 
perturbation. The details of these perturbations are all explained in detail in 
the text. The stationary state of these five cases correspond to a state in each 
metastable branch appearing in Fig. |3] although the branch corresponding to (a) 
does not appear in the numerical simulations with random initial conditions. 
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Figure 5. Velocity distributions in the case p = 2/7, corresponding to the states 
given in fig. |3] 

Next let us calculate the flow-density relation for each branch. In the metastable 
branches we find phase separation into a free-flow and a jamming region. In the former, 
pairs move with velocity t;/ and a headway of empty cells between consecutive pairs. 
In the jammed region, the velocity of the pairs is Vj and the headway dj. Nj and Nf 
are the numbers of cars in the jamming cluster and the free uniform flow, respectively. 
We assume Nf and Nj to be even so that there are Nf/2 and Nj/2 pairs, respectively. 
Then the total number of cars A'' is given by TV = Nj + Nf and the total length of 
the system becomes / = {dj + 2)Nj/2 -\- {df + 2)Nf/2. Since the average velocity is 
V — {NfVf + NjVj)/N and density and flow of the system are given by p = N/l and 
Q = pv, we obtain the flow-density relation as 

^^^-{^■-^■-'^- 

From the stationary states in Fig. 0] we have 

(a) : {vf,Vj,df,dj) = (5,4,6,4) 
(6) : {vf,Vj,df,dj) = (5,3,7,3) 

(c) : ivf,v,,df,dj) - (5,2,8,2) 

(d) : {vf,Vj,df,dj) = (5,1,9,1) 
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{e):{vf,v^,df,dj) = {5,0,10,0) (14) 

Therefore the resulting equations for each branch are 

Q = l + cp (15) 

where c= 1,1/2,0,— 1/2,— 1, which correspond to the branches A,B,C,D and E in 
Fig. El respectively. End points of the branches are found to be given by pi = 2/{df+2) 
and p2 = '2/{dj + 2), where pi is the point that the metastable branches intersect the 
free flow branch, and p2 is the maximal possible density in the metastable branches. 
Note that all p2 lie on the line Q = —2/3 + 2, which is indicated as the broken line in 
Fig. El 




Figure 6. A schematic diagram of the metastable branches (A, B, C and D) and 
the jamming Une {E) in the new model. The highest flow state is represented by 
Po, which is quite unstable and easy to go down to the lower flow state Pj\, ■ ■ ■ , 
according to the magnitude of the perturbation. 

Next let us now study the stability of each metastable branch. We mainly consider 
the density p — 2/7 and, in particular, wc will focus on the uniform flow represented 
by ■ • • 11000001100000 • ■ •, which shows the highest flow given in the point Pq in Fig. El 

Spatio-temporal patterns due to various kinds of perturbations are already seen 
in Fig. 01 Perturbation in this case means that some cars are shifted backwards at the 
initial configuration. The initial conditions for Fig.01[a)-(e) are given as follows: 

(a) Very weak perturbation (one car is shifted one site backwards)- ■ • 11000010100000 ■ • 

(b) Weak perturbation (one car is shifted two sites backwards) - - ■ 11000100100000 ■ - - 

(c) Moderate perturbation (one car is shifted three sites backwards)- - - 11001000100000 
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(d) Strong perturbation (one car is shifted five sites backwards) - • ■ fffOOOOOf 00000 ■ • ■ 

(e) Strongest perturbation (three cars are shifted backwards) ■ • • 11111000000000 • • ■ 

The stationary state of (a), (b), (c), (d) and (e) are given by the points Pa, Pb, 
Pc, Pd and Pe in Fig. IHl respectively. That is, if the system in Pq is perturbated, 
then the flow easily goes down to a lower branch in the course of time depending 
on the magnitude of the perturbation. Since the density does not change due to the 
perturbation, we obtain Pa ■ (2/7,9/7), Pb : (2/7,8/7), Pc ■ (2/7, 1), Pd : (2/7,6/7) 
and Pe : (2/7, 5/7) by substituting p = 2/7 into eq. lfT3|) . 

We see a jamming cluster propagating backwards in the cases of (d) and (e) in 
Fig. El In other cases the jamming cluster propagates forward ((a) and (b)) or does 
not move (c). These facts are related to the gradient of the metastable branches which 
are given by c according to eq. (|15|) . 

4. Stochastic generalization 

Finally we will combine the above model with the NS model in order to take into 
account the randomness of drivers. The NS model is written in Lagrange form as 

^x\+ max {0, min{F, x\+^ -x\-l,x\- x\-^ + 1} - r,*} . (16) 

where 77* = 1 with probability p and r/* = with probability 1 — p. The last term in 
the mininum in H16|l represents the acceleration of cars. The randomness in this model 
is considered as a kind of random braking effect, which is known to be responsible for 
spontaneous jam formation often observed in real trafSc j^j- We also consider random 
accerelation in this model which is not taken into account in the NS model. 

Thus a stochastic generalization of the hybrid model in the case of S* = 2 is 
similarly given by the following set of rules: 

1 Random accerelation 

=min{-l/„a.,«r^ +?7a}. (17) 
where rja — 1 with the probability pa and ?7a = with 1 — pa- 

2 Slow-to- accelerate effect 

(2) • r (1) t-i t-i CI 

3 Deceleration due to other vehicles 

(3) . I- (2) f t oT 

v\ = imn{v\ \ x,^s -x,- S\. 

4 Random braking 

(4) r (3) n1 

v] = max|Uj> - 77b,0) 
where r/f, = 1 with the probability pf, and rn, — Q with 1 — Pb 

5 Avoidance of collision 

.f+^)=min{.l"\x*+i-x*-l + (21) 

with n > 4, which is an iterative equation that has to be applied until v converges 

. (n+i) (") \ 

to v} = vl '(= Wj). 



(18) 
(19) 
(20) 
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6 Vehicle movement 



Vi 

„(0) 



(22) 



Again the velocity Vi is used as vl in the next time step. Step 5 must be applied to 
each car iteratively until its velocity does not change any more, which ensures that 
this model is free from collisions. This is the difference between the deterministic and 
stochastic case. In the deterministic model it is sufficient to apply the avoidance of 
collision stage only once in each update, while in the stochastic case generically it has 
to be applied a few times in order to avoid collisions between successive cars. 
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Figure 7. Fundamental diagrams and typical spatio-temporal patterns of the 
new stochastic model with different value of random parameters. Parameters are 
set to Vmax = 5 and 5 = 2. Upper two figures are the case of pa = and 
Pi, = 0.2, middle ones are pa = 0.8 and pi, = 0, and the bottom ones are pa = 0.8 
and pf, = 0.2. 

The fundamental diagrams of this stochastic model for some values of Pa and 
Pi, are given in Fig. \7\ The randomization effect can be considered as a sort of 
perturbation to the deterministic model. Hence some unstable branches seen in the 
deterministic case disappear in the stochastic case, especially if we consider the random 
braking effect as seen in Fig. [3 Random accerelation itself does not significantly 
destroy the metastable branches. Moreover, from the spatio-temporal pattern it is 
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found that spontaneous jam formation is observed only if we allow random braking. 
Random accerelation alone is not sufficient to produce spontaneous jamming. We also 
note that a wider scattering area appears if we introduce both random accerelation 
and braking. 

5. Concluding discussions 

In this paper we have proposed a new hybrid model of traffic flow of Lagrange type 
which is a combination of the BCA and the s2s model. Its stochastic extension is also 
proposed by further incorporating stochastic elements of the NS model and random 
accerelation. The model shows several metastable branches around the critical density 
in its fundamental diagram. The upper branches are unstable and will decrease its flow 
under perturbations. It is shown that the magnitude of a perturbation determines the 
final value of flow in the stationary state. Moreover, introduction of stochasticity in 
the model makes the metastable branches dilute and hence produces a wide scattering 
area in the fundamental diagram. We would like to point out that this metastable 
region around the phase transition density is similar to so-called synchronized flow 
proposed by j25f . Our investigation shows that one possible origin of such a region 
is the occurance of many intermediate congested states near the critical density. If 
some of them are unstable due to perturbation or randomness, then a dense scattering 
area near the critical density is formed around the metastable branches. This is in 
some sense in between the two cases of a fundamental diagram based approach (with 
unique flow-density relation) and the so-called 3-phase model of [2^1 which exhibits a 
full two-dimensional region of allowed states even in the deterministic limit. 
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